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Abstract 

A compact approximate groundstate of the Kondo problem is in- 
troduced. It consists of four Slater states. The spin up and down 
states of the localized d-impurity are paired with two localized s- 
electron states of opposite spin. All the remaining s-electron states 
are rearranged forming two new optimal orthonormal bases. Through 
a rotation in Hilbert space the two localized states (and the rest of 
the bases) are optimized by minimizing the energy expectation value. 
The ground-state energ y £"00 and the singlet-triplet excitation en- 
ergy AEst are calculated numerically. Although the two energies can 
differ by a factor of 1000, they are obtained simultaneously. The 
singlet-triplet excitation energy AEst is proportional to exp [—l/2Jp] 
and quite close to the Kondo temperature ksTx- The cases for anti- 
ferromagnetic (J > 0) and ferromagnetic (J < 0) coupling are inves- 
tigated. 

PACS: 75.20.Hr, 72.15.Rn 



1 Introduction 

The properties of magnetic impurities in a metal is one of the most intensively 
studied problems in solid state physics. The work of Friedel [T] and Anderson 
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[2] laid the foundation to understand why some transition-metal impurities 
form a local magnetic moment while others don't. Kondo [3] showed that 
multiple scattering of conduction electrons by a magnetic impurity yields 
a divergent contribution to the resistance in perturbation theory. In the 
following three decades a large number of sophisticated methods were applied 
to better understand and solve the Kondo model, and it was shown that at 
zero temperature a Kondo impurity is in a non-magnetic state. To name 
a few of these methods: scaling [1], renormalization [5], [6], Fermi-liquid 
theory [Zj, [8], slave-bosons (see for example [9]), and the Bethe-ansatz [10], 
[TT] . For a review see [12]. For numerical calculation an approximate wave 
function by Varma and Yafet [13] which was later extended to the large-spin 
limit [2], [15] was particularly productive. Finally after decades of research 
exact solutions of the Kondo and Friedel-Anderson problems were derived 
ilOj, il6j. 

One of the authors [T7], [18] recently introduced a new approach to the 
Friedel-Anderson impurity. This approach is based on building from the 
s-electron spectrum a localized s-state for each spin which behaves as an 
artificial Friedel resonance (AFR) states. This permits the construction of a 
very compact approximate ground state (see appendix, equ. ([2])). It requires 
solely the optimization of the two localized AFR states and consists of a 
few Slater states (Slater determinants). This solution gives remarkably good 
results for the ground-state energies and the occupation of the d-states. The 
results are of the same quality as numerical calculation by Gunnarsson and 
Schoenhammer p3] who applied the large Nf method to the spin 1/2 Friedel- 
Anderson model and included self-consistently up to 10'' basis states. 

In this paper we introduce a similar ansatz to the Kondo effect. The 
intension is to rederive some of the well known properties of the Kondo effect 
within a rather simple and transparent frame work. The hope is to extend 
this rather simple method to more complex problems such as the interaction 
between Kondo impurities. 

2 The AFR-Groundstate for the Kondo Ef- 
fect 

The Kondo ground state which is suggested in this paper is derived from the 
Friedel resonance problem. Below equation ([2]) gives the Friedel-Anderson 
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Hamiltonian. If one removes the Coulomb exchange interaction term Uud+nd- 
then one obtains the Hamiltonian for the Friedel resonance. This is a single- 
particle Hamiltonian and the n-electron ground state can be written as a 
simple product of single electron states which are a hybridization of the d- 
electrons with the band s-electrons. 

Due to the s-d-interaction between the s-electrons and the d-impurity the 
s-electrons shift their weight towards the d-impurity. This is more clearly 
expressed in a different version of the n-electron ground state of the Friedel 
impurity [20] . The system forms from the s-band a state which is localized 
in the close vicinity of the (i*-state. This Og-state hybridizes with the d*- 
state in the form {Auq + Bd*), where A,B are coefficients. The remaining 
s-electrons are arranged in states a*(l<i<A^ — 1) which are orthogonal 
to and each other. (Furthermore the matrix elements of the free electron 
Hamiltonian Hq must be diagonal (a*$o |-f^o| Oj^o) = E {i) 5ij for z,j > 1, 
see appendix). Then the alternative version of the {n + l)-electron Friedel 
ground state is given by 

i^Fr = {Aal + Bd*) \0a) (1) 

where the hybrid {Aa^ + Bd*) plus the n lowest states 0^03. .a* = \0a) of the 
basis {a*} are occupied. The state Oq can be found by numerical variation 
or from an analytic expression (from Oq the full basis {a*} can be derived). 
One of the authors has proved [2T] that ([T]) is an exact ground state of the 
Friedel impurity. 

In the next step we consider the Friedel- Anderson impurity whose Hamil- 
tonian is given in equ. ([2]). 

N N 

Hfa = '^{'^eucl^Cua+Eddld^+'^Vsd{i^)[dlcua+cladcT]}+Und+nd- (2) 

CT U=l U=l 

Its mean field solution is just a product of two Friedel solutions as given by 
equ. ([1]), one for the spin up and the other for the spin down electrons (with 
different effective d*-state energies Edi and Edi). It requires two localized 
s-states, which we denote now as for spin up and 6q for spin down and 
the rest of their bases {a*} and {b*} . The two differ when the impurity is 
magnetic. 

The mean field solution can be dramatically improved if one expands 
the product of the two hybridized states [A-^a^^ + B^d'^) {^iKi + ^l^l) 
assigns new coefficients {A, B, C, D) for four two-electron states. This yields 
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ijms = [Aa*^bl^ + Ba*^dl + Cd*bl^ + Dd*dl\ \0a^) \Obi) (3) 

The optimal state ipms is obtained by varying the composition of Oq and 
6q and optimizing the coefficients {A, B,C, D). This solution has the same 
structure as the mean field solution but it yields a much lower ground-state 
energy. It also requires a much larger threshold Coulomb exchange energy U 
to form a magnetic moment at the impurity. 

The study of the Kondo effect taught us that the real ground state of 
the Friedel-Anderson impurity is a symmetric state (often called a singlet 
state). Such a state can be constructed from ([3]) by reversing the spins of all 
electrons and combining the two states. This yields 

^lJ,s = [Aa*^b*^ + Bal^d\ + Cd]hl^ + Dd^dj] |0„t) (4) 
+ [A'bl^al^ + C'bl^dl + Bd]'al^ + D'd]d*^ |Obt) 10^;) 

In equ. (|4]) the creation operators have been reordered according to their 
spin. By optimizing the localized states and the coefficients we obtained 
an almost perfect agreement with Gunnarsson and Schoenhammer (TU] for 
the ground-state energy and the zero, single and double occupation of the 
(i*-state. 

In this paper we are interested in the Kondo effect. For the treatment 
of a magnetic impurity Kondo used the exchange Hamiltonian Hsd with the 
exchange interaction J^^k'- Oiig generally approximates the exchange inter- 
action by a 5 (r)-function: J (r) = Va8 (r) where Va is atomic volume. Then 
the exchange (or Kondo) Hamiltonian has the form 

r (^+vI;t(0)vI/^(0) + 5_vI/t(0)vI;^(0)) 
[ +S. (vl/| (0) (0) - vl/t (0) vp^ (0)] 

where S+, S-, Sz are the spin operators of the impurity with spin 5 = 1/2 and 
(0) and (0) represent field operators. The product VaJ"^^ (0) \I'o-' (0) 
yields an energy since (0) (0) has the dimension of a density. 

Schrieffer and Wolff [22| showed that there is an intimate connection be- 
tween the Friedel-Anderson and the Kondo impurity. The Friedel-Anderson 
Hamiltonian can be transformed into an exchange Hamiltonian with an anti- 
ferromagnetic exchange interaction J > (plus some additional terms). In 
particular for large values of U and large negative (for example E^, = —U) 



(5) 
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they derived an effective J 



The Kondo impurity is in some respects a hmiting case of the Friedel- 
Anderson impurity. By increasing the exchange interaction U (and decreas- 
ing Efi = —U/2) the Friedel- Anderson impurity approaches the properties 
of a Kondo impurity. For the Kondo impurity the d-impurity is a locahzed 
permanent magnetic moment. This means that the (i*-state is always singly 
occupied, either with spin up or down. We observe this transition in our so- 
lution ipss (equ. H]). When we increase U (with Ed = —U/2) the coefficients 
A, A' and D, D' approach zero. It becomes energetically to expensive to have 
an empty or doubly occupied (i*-state. Therefore for the Kondo impurity we 
make the following ansatz for the ground state 

^l;K = {Ba*^dl + Cd}b*^) |0,^0,i) + {Cb*^dl + Bd*a*^) \0,^Oai) (7) 

= B'ipB + Ci)c + Cipc + 

We remove the terms with zero and double d-occupancy from the AF ground 
state ipss- 

We denote the state ipx as the AFR ground state of the Kondo effect 
since it uses the artificial Friedel resonance states and 6q. 

As in the AF ground state ipss the states Oq and 6q are localized s-states 
which are concentrated close to the impurity. Their composition will be nu- 
merically optimized. Here, we sketch briefly the construction of the {aQ,a*} 
basis. (The construction of the {6q, h*} is completely analogous). 

In this paper we use a finite s-electron band with states c* per spin 
and a magnetic d*-sisXe with spin S = 1/2. For the present calculations we 
use Wilson's band which extends from —1 to +1 and has an exponentially 
fine energy spectrum close to the Fermi energy (see appendix). The ag-state 
is composed of these c*-states 

N 

«o = 5Z (8) 

u=l 

(In the initial phase the coefficients are arbitrarily chosen, for example 
being constant). 

The remaining s-electron basis{c*} has to be made orthogonal to a^. 
This yields the new remaining bases {a*}- The a* {1 <i < N — 1) are built 
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orthogonally to and to each other. In addition their (A^ — 1) x (A^ — 1) 
sub-matrix of the s-band Hamiltonian Hq = ^u^u is made diagonal. The 
a* together with represent a new basis. (Details of the construction of the 
states {oq, a*} is discussed in ref. [17], [18] and briefly in the appendix). The 
states a* are uniquely determined from the state Oq. Their form is 

N 

The new bases {a*Q,a\) can be expressed as a rotation of the original basis 
{c*} in Hilbert space. The basis {feo;^!} is built in complete analogy. 

The free electron Hamiltonian can be expressed in either of the two new 
bases. It is diagonal in the states a* for 1< z < (A^ — 1) and has the form 

Ho = ^ Evclcy = ^ E (?) alai + E (0) a^a^ + ^ Vf^ {%) [a^ai + a* oq] (9) 

v=\ 1=1 i=l 

The Hamiltonian is equivalent to a Friedel Hamiltonian where a*^ is an 
artificial Friedel resonance state (APR state) which is purely composed of 
s- states. 

After the initial construction of the two bases {aQ,a*} and {&05^i} the 
APR ground state ([7]) of the Kondo impurity can be formed. 

3 Numerical Results 

The procedure to obtain the optimal states Oq and feg has been described 
in previous papers [20], [17], [H] and is briefly sketched in the appendix. 
One starts from an s-band with A^ states possessing the energy Sy. In all 
calculations the energy band ranges from —1 to 1 following Wilson's example. 
The density of states is constant and equal to p = 1/2 (corresponding to 
one state in the whole band). This band is divided into A^ adjacent cells 
{Ey-i : Ejj). The energy values in the center (middle) of each cell yield the 
spectrum = {Ey — E^^i) /2. 

It was pointed out by Wilson and will be confirmed in the present calcu- 
lation, that a very small energy spacing 6E at the Permi energy is of essential 
importance for obtaining the Kondo ground state. A linear energy scale is 
therefore not well suited for the Kondo ground state because SE is always 
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much too large. We use an energy spectrum with an exponential energy 
scale. This energy spectrum was introduced by Wilson [5] in his Kondo pa- 
per. Starting with the values of —1 and +1 the energy is N/2 times divided 
by 2, yielding the cells (-1 : -1/2) , (-1/2 : -1/4) , .. (-1/2^/^ . g). The 
negative value of are (for I < v < N/2): -3/4, -3/8, .., -3/2''+\.. - 
2^2^/2+1^ —1/2^/^+^. The positive energies are the mirror image of the neg- 
ative ones. 

From the chosen energy spectrum the two states Oq and h*Q are constructed 
with initially arbitrary coefficients and /3q. The states and deter- 
mine uniquely the full bases {a^^a*} and {hl^h*} (see appendix). With the 
two bases the energy matrix elements between the different Slater states are 
calculated. 

The wave function has four components and the Hamiltonian consists 
of four terms due to Hgd and two terms due to Hq. Therefore the energy 
expectation value is composed of 96 matrix elements. Many of these are 
identical and one has to calculate onlyl6 different matrix elements. They 
can be transformed into determinants of (ra + 1) x (n + 1) matrices. In the 
appendix some representetive matrix elements are derived. 

Since the two Slater states {ipB, V'c) ^'^^ orthogonal (nor are {ipc, V^)) 
one has first to apply an transformation to the four-component basis (5, C, C, i?) 
before one can solve the eigenvalue problem. (Details of the optimization are 
described in the appendix.) 

The lowest value of A yields the lowest energy expectation value {E) for 
the chosen two bases {oq, a*} and {6q, h*}. Then states and 6q are rotated 
in Hilbert space until the value {E) reaches a minimum. The resulting state 
is defined as the AFR ground state of the Kondo impurity and the energy 
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Eqq = (E) is its ground-state energy. 




Fig.l: The energy of the AFR ground state of the Kondo 
impurity as a function of the exchange interaction J for 
= 40 states. 

This ground-state energy Eqq is plotted in Fig.l as a function of the 
exchange interaction J. The curve shows the results for the exponential 
spectra with = 40. The corresponding value for the energy spacing at the 
Fermi energy is ^-E = 1/2^° » 1 x 10^^. 

The ground-state energy Eqq can be well expressed by the following de- 
pendence 

Eqo (J) = -.5017 * J2 + .435 * - 9.65 * jH 

This ground-state energy is not the Kondo energy. It is sometimes called 
the perturbational part of the ground-state energy. It is (in most cases) 
ignored in the Kondo effect as the non-singular part of the ground-state 
energy. Below we will compare this ground state with the ground-state energy 
for ferromagnetic coupling, i.e. for negative values of J. 

The solution of the eigenvalue equation f|T^ yields four eigenvalues. The 
lowest eigenvalue corresponds to the ground-state energy. This should be 
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a singlet state. We expect that the next eigenvalue of the energy, corre- 
sponding to the first excited state, represents a triplet state. To confirm 
these expectance we calculate the expectation value of the total spin squared 

= (^jSj)^ for the two states. For J = 0.1 we find in the ground state 
(S^) — 0.04. For the first excited state we obtain for the square of the total 
spin (S^) = 1.99. This is very close to the value 2 which we expect for S — 1. 
We conclude that the ground state is not a perfect singlet state but is pretty 
close to it. It is essentially a state with = 0, but it has a small component 
of a larger total spin such as S = 1,2... Similarly the first excited state has 
essentially the total spin S — 1. But it too has a small admixture of states 
with S equal to 0, 2, ... (The fact that (S^) is so close to 2 means that the 
contributions of S* = and S = 2 essentially cancel each other). For the rest 
of the paper we denote the ground state as the singlet state ipK,o- The first 
excited state we denote as a triplet state ipK,i- 

The energy difference between these two lowest energy states we denote as 
the singlet-triplet excitation energy AEgt- This energy difference A.Est is very 
small and In (AEgt) is shown in Fig.2 in a logarithmic plot as a function of 
the inverse coupling strength 1/J. AEgt follows an exponential law covering 
six orders of magnitude in the range 0.05 < J < 0.18. The linear full curve 
in Fig.2 is given by AE w 5 * D * exp [—x] where x = 1/ (2Jpo)- The other 
data and curves in Fig.2 are explained in the discussion. 
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Fig.2: The energy difference between tfie singlet and triplet 
state. The full circles are obtained as the energy difference 
between the ground state and the first excitated state. The 
stars and the dashed and dotted curves are explained in the 
discussion. A exponential Wilson band with A'^ = 40 
and A'" = 60 is used in the calculation. 

Wilson has shown that the energy separation of the states at the Fermi 
energy has to be smaller than the Kondo energy to obtain the full Kondo 
effect (or in Wilson's words: to make the transition into the infinitely strong 
coupling case). This same behavior is observed in our APR solution. 

Wc investigate how critical the smallest level separation 5E at the Fermi 
energy is. For this purpose we vary the number A^ of states for the expo- 
nential energy spectrum. The level spacing at the Fermi energy 5E depends 
exponentially on A^: (5£'=2^^/^.For each A^ we calculate the AFR ground 
state for two values of J, J = 0.1 and J = 0.07. In Fig. 3 the singlet-triplet 
energy difference ^Egt is plotted versus the smallest energy spacing 5E in a 
log-log plot. The number A^ of energy levels is varied from A^ = 12 (at the 
right side) to A" = 48 (left side). One recognizes that for sufficiently small 5E 
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(on the left side of the figure) the singlet-triplet excitation energy assumes 
a constant value. With increasing 5E the excitation energy increases. The 
dashed straight line has a slope of 0.85 corresponding to an increase of /S.Est 
proportional to {6E)^'^^ . Fig. 3 demonstrates that the transition to a constant 
AE'st occurs roughly for 5E w ^Egt and that 5E must be at least a factor 
10 smaller than the final result for AE'^t to give a rehable value for l^Egt- 




— I 1 — ■ 1 ' 1 ' n 

-16 -12 -8 -4 

ln(5E) 



Fig.3: The logarithm of the singlet-triplet excitation energy versus 
the logarithm of smallest energy spacing 5E at the Fermi energy. 



3.1 AFR ground state for ferromagnetic coupling 

One of the interesting aspects of the exchange Hamiltonian is that it only 
yields the Kondo anomaly for anti-ferromagnetic coupling or J > 0. We 
apply the AFR-ansatz to the case of ferromagnetic coupling, i.e. J < 0. In 
Fig.4 the resulting ground-state energy is plotted versus the absolute value 
|J|. The ground-state energy lies somewhat above that of the Kondo case 
with anti-ferromagnetic coupling. The difference is about 10% for small \J\ 
and a factor two for \J\ w 0.18). 
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It is in particular interesting that the two lowest eigenvalues of equa- 
tion (1141) differ on by about 10~^^ for sufficiently small level separation SE. 
Here the value of |Ai?st| is always much smaller than 6E. In other words, 
at all temperatures the two states, singlet and triplet, are degenerate and 
the impurity is magnetic. The AFR solutions for anti-ferromagnetic and 
ferromagnetic coupling are distinctively different. 




Fig.4: The ground-state energy of the system for negative J, 
i.e., for ferromagnetic coupling. 

Eo^f = -.4630 * J2 + .680 * - .697 * 

3.2 The artificial Friedel resonance states 
3.2.1 Anti-ferromagnetic coupling (J > 0) 

The whole bases {aQ,a*} and {bQ,b*} can be derived from the composition 
of the AFR states Oq and b^. In Fig.5a,b the coefficients and Pq of the 
states Oq and feg plotted versus the state index u ior N = 40 (using 
equ.®). On the left and right side of u = N/2 the values of u and (A^ — u) 
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represent essentially the logarithm of the energy. The coefficients are always 
positive because this yields the maximal interaction of the AFR states with 
the magnetic d* state. The state Oq is mainly composed of states below sp 
and has its main weight at energies above Sp- However, in the center close 
to the Fermi energy the amplitudes of the two locahzed states are essentially 
identical. This is a behavior which we also observed for the Friedel- Anderson 
impurity. 
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Fig.5a,b: The coefficients of the AFR states ckq and /Sq as a function of i/ for the 
two AFR states Oq and h^. The exchange interaction is for anti-ferromagnetically 
coupling (J = 0.1). 



For large energies the amplitudes and (3^ of Oq and 6q are rather dif- 
ferent. For the analysis at small energies we plot the occupation density 
Ictgl^ / {Ej, — E'lz-i) and l/Sg |^ / {Ey — E'v-i) as a function of u. This is shown 
in Fig. 6 where both occupations are plotted in the same figure. At energies 
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close to the Fermi energy the occupations of Oq and 6q are almost identical. 



800- 




Fig.6: The occupation of the and APR states for = 60. Note that the 
figure shows two curves as a function of u {0 < u < 60). At small energies 
(close to the center arrow) the two are almost identical. 

3.2.2 Ferromagnetic coupling (J < 0) 

For comparison we plot the coefficients and Pq for J = —0.1, i.e. the 
case of ferromagnetic coupling in Fig. 7a, b. For the ferromagnetic coupling 
we do not observe the two maxima as in the Kondo case. Furthermore, 
the coefficients have negative sign in half of the energy range, reducing the 
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amplitudes of the AFR states at r = 0, the posit 



ion of the (i*-moment. 
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Fig.7a,b: The coefficients of the AFR states and /3q as a function of u for the two 
AFR states and b^. The exchange interaction is ferromagnetically couphng (J = —0.1). 



4 Discussion 

Yosida [23] introduced in the 1960's an approximate solution for the Kondo 
Hamiltonian of the form 

^ = ^(«OT^I + «i)|0) 

where 

io) = n 4,.$o 

cr,k 

k<kp 

is the quasi-vacuum with all states below the Fermi energy occupied and 
Oq is an optimized single electron state which is composed of states with 
k > kp. The state Oq is therefore orthogonal to all occupied states. $o is the 
vacuum state. Yoshida's ground-state energy and singlet-triplet energy were 
of the same order of magnitude with 

£"00 = D exp 



fJp 
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Yosida's ansatz and our suggested AFR ground state are related to each 
other in may respects: (i) Both use a locahzed state Oq, (ii) both form pair 
states between the locahzed state and the d* state with zero spin component 
in z-direction, (iii) both fill the remaining electrons into states which are 
orthogonal to Cq. The main difference is that in our approach the state 
together with a second state are composed of all band states. This requires 
the remaining electron states to be aligned perpendicular to the localized 
states. Our AFR ground state is in a way a revival of Yosida's approach. 

4.0.3 Ferromagnetic versus anti-ferromagnetic coupling 

The numerical calculations show that the ground state for the anti-ferromagnetic 
(J > 0, Kondo effect) and the ferromagnetic coupling (J < 0) are very dif- 
ferent. 

• The groimd-state energy for positive J lies lower. The behavior is 
roughly described by the two fits to the numerical results: 

- For J > 0: Eq+ = -.5017 * + .435 * - 9.65 * 

- For J < 0: E^- = -.4630 * | + .680 * | - .697 * | 

• There is essentially no singlet-triplet excitation energy in the case of 
ferromagnetic coupling. The reason is that the system adjusts in such 
a way that not only the states with opposite d*-state spin (for example 
■0B and ■0c) are orthogonal but even the states with the same d*-state 
spin orientation (for example iI)b and ip-^) are quasi-orthogonal. This 
becomes visible in the (multi-electron) scalar products between the 
state iPb,4'ct'4'c,'4'b- ^ — ~0.1 the scalar products {iPb\'4'c) 
{ipc\i^) are of the order of 10~^ while in the anti- ferromagnetic Kondo 
(J = 0.1) {'4>b\'4^) and are of the order of 1/3. In the case 
of the ferromagnetic coupling J < the two bases (oq, a*} and {6q, 6*} 
align themselves in such a way that the scalar product between |0a) 
and |0a) with n occupied states is of the order of (0a|0b) ~ 10""^. The 
bases align themselves quasi-orthogonal. Therefore the reduction of 
the ground-state energy is solely due to z-component of the interaction. 
The energy matrix element of the x-y-component of the interaction is 
practically zero. 
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4.0.4 Coefficients of the ground- and excited state 

Our Kondo-ground state ipK is given by ([7j). Tlie coefficient- vector {V) = 
(5 C C B y contains interesting information about the structure of 
the states. Here we discuss the resuhs for J = 0.1, = 40. In the ground 
state {Vy takes the values (^)^,ound = (•703,.0626,.0626,.703) and in the 
excited state we find iV)l^^-^^^ = (.705,-0635,0635,-705). (Note that state 
ipK,o and ipK,i are normalized but not the vectors {V) because the pairs of 
states ipB,ip(j and ipci^ are not orthogonal as discussed in the appendix). 
For both, ground and excited states, the sub-states with Oq have a large 
amplitude while the sub-states with the state 6q have a small amplitude. In 
the ground state all coefficients are positive. In the excited state the sub- 
states with opposite (i*-spin have opposite sign. These are the signatures of 
a singlet ground state and triplet excited state. (Note that the states are 
ordered according to their spin). 

For the ferromagnetic coupling we obtain for J = —0.1, = 40 the 
following coefficients: In the ground state {V^ takes the values (V")^^,^^^^^ = 
(-.706, .-.0400, .0400, .706). For the essentially degenrated excited state we 
find iV)l^^-^^^ = (.706,.-.0400,-.0400,.706). This means that in the ground 
state the components with BaQ^d* and Bd'^aQ^ form a triplet state (as do the 
components C(i|6Q| and Cbl-^d*^) but the signs of the states which contain 
and 60 with the same spin are opposite. 



4.0.5 Singlet-triplet excitation energy 

The triplet state ipK,i is obtained by minimizing the ground-state energy. 
Then ipK,i is the first excited state. This state ipK,i is composed of the same 
Slater states ipB,4'c,4'~B,4'c singlet state. Only the coefficients, and 

particularly their signs, are different than in ipK,o- But this state is not 
necessarily the triplet state with the lowest energy. To obtain the lowest 
triplet state during the minimization we impose on the coefficients B, C and 
B, C the conditions 

B = -B C = -C 

This automatically enforces the triplet state. With these conditions we re- 
peat the optimization procedure to obtain the lowest triplet energy. The 
resulting state is a relaxed triplet state. We denote it as ipK,i- Its Slater 
states ipB, i'Ci i^, i^c different than for the singlet ground state i^xfi- Its 
energy Ek,i is lower than the energy Ek,i of the unrelaxed triplet state ipK,!- 
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Now we can take the energy difference between the energy of the relaxed 
triplet state and the singlet state. Because the two values are obtained in two 
different optimization calculations the absolute accuracy of the optimization 
has to be better than 10~^° (in units of the band width). We plot the new 
excitation energy in Fig.2 as stars. It lies below the unrelaxed excitation 
energy. The latter is approximately given by AE^t ^ 5Dexp[— 1/ (2Jpo)]- 
Our new values for the excitation energy lie between two theoretical curves: 
(i) AEst = -Dexp [—1/ (2Jpo)] , given by the dashed curve and (ii) AEgt = 
v^2JpoZ}exp [—1/ (2Jpo)], given by the dotted curve. Both expressions are 
given in the literature as approximate values for the Kondo temperature 
ksTK- The numerical values lie closer to the second expression. 

4.1 The artificial Friedel resonance states 

In the Fig.5a,b and 7a,b the coefficients (or /3q respectively) as they occur 
in equ. ([H]) are plotted as a function of u. One has to recall that for = 40 
the Fermi energy lies in the center at z/ = 20. The energies above the Fermi 
energy are: i2~^^, |2~^^, |2~^^, ..|2~^ The energy states below the Fermi 
energy are a mirror image of the state above. 

The Kondo and the ferromagnetic case have in common that the coef- 
ficients of and fcg ^^e almost mirror images. The state Oq has its main 
weight at negative energies and at positive energies. 

Beyond that the ferromagnetic and anti-ferromagnetic cases differ con- 
siderably. The coefficients in the Kondo case are all positive whereas in the 
ferromagnetic case the signs in the positive and negative energy regimes are 
different. 

5 Conclusion 

In this paper we have developed a compact (approximate) ground state for 
the Kondo problem. Two localized s-electron states, and 6q are built from 
the s-band. Their amplitudes, normally two times 40 numbers, determine 
fully the ground state. This ground state consists of four Slater determinants. 
It yields two energies, the total ground-state energy Eqq and the singlet-triplet 
excitation energy AEst- While Eqq depends essentially quadratically on the 
coupling constant J, the energy AEst shows an exponential dependence. It 
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is given by 

and is proportional to the Kondo temperature. 

The energy scales of AEst and Eqq are very far apart. They differ by a 
factor 100 to 1000. It is remarkable that our solution works for both energy 
ranges at the same time. 

Our ansatz shows also a very different solution for the Kondo (J > 0) 
case than for the ferromagnetically coupled case (J < 0). In particular the 
latter has no singlet-triplet splitting, and the different sub-states ipx avoid 
each other by forming quasi-orthogonal multi-electron states. 

The structures of the locahzed AFR states and 6q differ considerably 
between the J > and the J < case. In particular for the J > case 
the coefficients have two maxima in an intermediate positive and negative 
energy range. A better understanding of these structures is left for future 
investigations. 

Acknowledgement: The research was supported by the National Science 
Foundation DMR-0439810. 
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A Appendix 



A.l Wilson's s-electron basis 

Wilson [Sj in his Kondo paper considered an s-band with energy values rang- 
ing from —1 to 1. In the next step Wilson replaced the continuum of s-states 
by a discrete set of states. This is done on a logarithmic scale. The discrete 
energy values are 1, 1/A, 1/A^, etc and —1, — 1/A, — etc where A is 
a parameter larger than one. (In this paper A = 2 is chosen). These dis- 
crete points are used to define a sequence of intervals: the interval z/ (for 
j/<N/2) is = -1/2"'^ < e < -1/2" = (there are equivalent inter- 
vals for positive ^-values where u is replaced by (A^ — z/) but we discuss here 
only the negative energies). The new Wilson states c* are a superposition 
of all states in the energy interval {^u~i,(.u) and have an (averaged) energy 

i^u ~\~ ^u-i) = ^~2^ 2^' ~4'~i'~^'"'~ 2-2"^/2 ' ~ 2-2'^/2 ■ This spec- 
trum continues symmetrically for positive energies. The essential advantage 
of the Wilson basis is that it has an arbitrarily fine energy spacing at the 
Fermi energy. 

While Wilson chose the logarithmic energy scale for his scaling approach, 
in this paper the physical background is somewhat different. We start (like 
Wilson) with the somewhat artificial energy band in the range (—1,1) which 
has a constant density of states. The volume of the host is Vh and the 
atomic volume va- The full band (per spin) shall have one electron per atomic 
volume. So the total number of electrons of a full band in the volume Vh is 
Z = Vh/vA- The electron density per spin of a full band is n = Z/Vh = 1/va- 
The density of states (per spin) is Z/2 in the volume Vh and equal to p = 1/2 
in the volume va- 

In this model the original energy levels have a (constant) separation of 
SE = 2vA/Vh which will be in the range between 10^*^ to 10~^^. In an energy 
interval {C,u~i,C,u) one has a large number Z^, of electron states f^{r) with 
1 < H < Zi,. Their level separation shall be constant and the value of the 
wave function at the origin ip^ (0) = l/y/Vi is independent of /i. 

Next we form a new basis Xa{r) for the energy range (^i., .^i^-i). We 
renumber all the original states y?^ (r) in this energy range from 1 to Z^, and 
define the new state (or wave function) Xa (r) (1 < « < Z^) 

Xa (r) = -^Et=i^, (r) e^^™^/^" 
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In the new basis Xq. ^ states except xz^, (^) have a zero amphtude at the 
origin. 

r =0 for \<a<Z^ 



for 



In the following wc denote the state xz^ i^u s-nd neglect all other Xa (r)- 
states because the latter do not interact with the impurity. The field operator 
^ (r) in the original basis ip^ has the form (c<^^ is the annihilation operator 
for the state ^pi, (r)) 



* (r) = E;=i<^m (r) c, 



We can divide ^ (r) into one part which is non-zero at r = and a second 
part which vanishes at zero. For the interaction with the d-impurity only the 
first part contributes. This part is equal to 



At r = one has 



* (r) = (r) c. 

*(0) = *(0) = EV^, (0) c, 
I ^ fz: 1 




/Va V Z 



The width of the energy ranges can be chosen arbitrarily. We use Wilson's 
exponential energy spectrum. This resulting basis {c*} is the starting point 
for all calculations. 

Since J is the matrix element per atomic volume we have to use also the 
density of states per atomic volume. The latter is 

1 

^=2 



A. 2 Construction of the Basis a, 
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a* 



For the construction of the state Og and the rest of basis a* one starts with 
the s-band electrons {c*} which consist of states (for example Wilson's 
states). The (i*-state is ignored for the moment. 



• In step (1) one forms a normalized state Qq out of the s-states with: 
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N 

E 

i/=i 



* 



(10) 



The coefficients can be arbitrary at first. One reasonable choice is q;° = 

1/Va^ 

• In step (2) {N — 1) new basis states a* {I < i < N — 1) are formed 
which are normahzed and orthogonal to each other and to Qq. 

• In step (3) the s-band Hamiltonian Hq is constructed in this new basis. 
One puts the state Oq at the top so that its matrix elements are i?oi 
and Hio. 

• In step (4) the {N — l)-sub Hamiltonian which does not contain the 
state Oq is diagonalizcd. This transforms the rest of the basis {a*} into 
a new basis {oq, a*} (but keeps the state Oq unchanged). The resulting 
Hamilton matrix for the s-band then has the form 



m 

Vfr{2) 



Vfril) 

m 





\VfriN-l) 



Vfr{2) 



E{2) 




Vfr{N 






1) \ 



(11) 



E{N -I) I 



The creation operators of the new basis are given by the set {oq, a*} , 
(0 < i < — 1). The a* can be expressed in terms of the s-states; 
a* = '^^c*. The state Oq uniquely determines the other states 

a*. The state Qq is coupled through the matrix elements Vfr (i) to the 
states a*, which makes the state an artificial Friedel resonance. The 
matrix elements E (i) and V/r (i) are given as 



V 

In the final step (5) the state is rotated in the iV-dimensional Hilbert 
space. In each cycle the state is rotated in the (aQ,a*J plane by an 
angle ioi 1 < iq < N — 1. Each rotation by Oi^ yields a new oq* 



dQ* — Oq cos Oi^ + a* sin 0, 
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The rotation leaves the whole basis {a^, a*} orthonormal. Step (4), the 
diagonalization of the (A^ — l)-sub Hamiltonian, is now much quicker 
because the [N — l)-sub-hamiltonian is already diagonal with the ex- 
ception of the if)- row and the zo-column . For each rotation plane 
(aQ,a*^) the optimal with the lowest energy expectation value is 
determined. This cycle is repeated until one reaches the absolute min- 
imum of the energy expectation value. In the example of the Friedel 
resonance Hamiltonian this energy agrees numerically with an accuracy 
of 10~^^ with the exact ground-state energy of the Friedel Hamiltonian 
|20j . For the Kondo impurity the procedure is stopped when the ex- 
pectation value changes by less than 10~^° during a full cycle. 

A. 3 Matrix elements 

In equ.([5]) we expressed the Kondo Hamiltonian in terms of field operators. 
We can also express it in the basis {c*}. Here it has the form 



with 



S+cl^ (0) c,,T (0) + (0) c,,^ (0) 
+S. (0) c.,T (0) - 4^ (0) c.,; (0) 



Jyy = J\j = - Cu-l 

We will, however, use (a slighly modified version of) equ.Q. The non- 
spin-fiip part can be rewritten as 

Hnsf = VaJS, (-^t (0) ^{ (0) + (0) ^{ (0)) (12) 

where we anti-commuted the field operators. 

We consider first the example of the non-spin-flip part where the spin of 
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the d-impurity and the s-electrons are anti-parallel: 

(^B I^ttI = -^'-^ {al^dlOa^Ob^ \S,^^ (0) ^{ (0)| a^^rfjO^tOfe^ 

_ VaJ_ 
~ 2 

v„J 



a*^dlOa^O,^ *T (0) (0) a*^c^*0„T06; 



= ^ (0) aotOaTC^lOH 1 1 ^{ (0) a*^0„tdjOH 



i=0 



where is the amphtude of the state a* at r = 0. The summation is from 
< i < n. Similarly one obtains for the non-spin-flip part where the spin of 
the d-impurity and the s-electrons are parallel: 



VaJ 



i=l 



A'^ is the amplitude of the state b* at the origin. The latter summation is 
from 1 < i < n. 

As an example for the spin-flip matrices one obtains 



(V^ \H^^ I iPb) = VaJ (d;a*^ObTO„i S+^\ (0) (0) al^dlOa^O,^ 

= -VaJ (^1 (0) d*^a*^0h^0ai II II (0) a*^d*^0a^0hi 
= -VaJ (*| (0) Ob^a*^Oai II II a^rOaT*! (0) Ob;) 

= J (a*^0,i 1 1 M/| (0) 0,^) (v]/t (q) q,^ | | a*^Oa^ 

^-VaJ\{^H0)Ob\\a*oOa)f 

The final brackets () represent a multi-scalar product which is given by the 
determinant of the following (n -|- 1) x (n -|- 1) matrix. 

(*t(o)o,||;o„> 

{^*{0)H) (**(0)|at) (**(0)|<) 

(blH) (blK) (KK) 



(KK) 



(KK) 
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^0 



AS 



mo 






(KK) 


(KK) . 


■■ (KK 



A. 4 Energy optimization 

For a given set of bases {a^, a*} and {60, b*} the energy expectation value of 
the Kondo state is given by 



The numerator of the expectation value of the ground-state energy is 
given by 



{i^Kmi^K) = {B C C B) 



where {V) = { B C C B ^ . The denominator is given by 



/ Hbb 


Hbc 


Hbc 


Hbb 


\ 




( B\ 


HcB 


Hcc 


Hcc 


HcB 






C 


HcB 


Hcc 


Hcc 








c 




Hbc 


Hbc 


Hbb 


) 







(13) 



(i^Kl^K) ^ (V)^ (N) (V) 



where (iV) is the non-diagonal symmetric matrix of the scalar products 
(since the Slater states V'x are pair- wise not orthogonal) 



(7V) = 



/ 1 (^bIV'c) 

{^\^b) 1 



\ 






1 






1 



With an orthogonal transformation of (V) the matrix (N) can be transformed 

into a diagonal matrix (A^) = (NaSij). Replacing V^— > / \/ transforms 
the matrix (N) into the unity matrix. 

Now we can require that the new vector (V) is normalized so that 

B'"^ + C + C"^ + B 1=1 and we have only to vary the nu- 
merator. This yields the eigenvalue problem 



{H') (V) = A (V) 



(14) 
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which can be easily solved numerically. The original coefficients B C C B ) 
are obtained by reversing the transformations. 
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